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In this paper we tackle the issue of causahty in quantum gravity, in the context of 3d spin foam 
models. We identify the correct procedure for implementing the causality /orientation dependence 
restriction that reduces the path integral for BF theory to that of quantum gravity in first order 
form. We construct explicitly the resulting causal spin foam model. We then add matter degrees of 
freedom to it and construct a causal spin foam model for 3d quantum gravity coupled to matter fields. 
Finally, we show that the corresponding spin foam amplitudes admit a natural approximation as the 
Feynman amplitudes of a non-commutative quantum field theory, with the appropriate Feynman 
\^ • propagators weighting the lines of propagation, and that this effective field theory reduces to usual 

' QFT in flat space in the no-gravity limit. 
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, In recent years, spin foam models Q, ^ |^ have attracted much interest as a new implementation of the 
ff^ ' sum-over-histories quantization of gravity |j| in a discrete setting and in a purely group-theoretic and algebraic 
' language. Spin foam models can in fact be understood as background independent discrete quantum gravity path 
integrals: spacetime (and its topology) is represented by a combinatorial 2-complex (collection of vertices, links 
I and faces) and a possible configuration of geometric data is determined by a collection of either (Lorentz) group 
elements or group representations assigned to its links or faces, respectively; each geometric configuration is assigned 
a quantum amplitude (usually, but not necessarily, derived via a discretization of some continuum formulation of 
y—i . classical gravity), and the partition function of the model is defined by the typical sum over all possible quantum 
OO ' configurations weighted by the chosen amplitude. The 2-complex is usually chosen to be topologically dual to a 
' simplicial complex of the appropriate dimension. The situation as regards 4-dimensional quantum gravity in a spin 
foam context is still in rapid evolution, with various models being proposed and several issues still to be clarified 0,12] ■ 
] On the other hand, the viability of a spin foam quantization of 3d quantum gravity is not anymore an issue, as it has 
O . been shown that the Ponzano-Regge spin foam model (see for example [3|) provides, after suitable regularization, a 
^ '' well-defined and useful description of it, whose equivalence with other approaches, e.g. Chern-Simons quantization 
$H I |13| or canonical loop quantum gravity 0, can be established. Moreover, spin foam models can be obtained from 
the perturbative expansion of 'group field theories'!^ 0, II], peculiar quantum field theories over group manifolds 
providing a kind of 'local simplicial third quantization of gravity', in the same way as the traditional path integral 
for quantum gravity can be obtained (formally) from the perturbative expansion of a suitably defined field theory 
on superspace (space of geometries). In doing so, one obtains also a definite prescription for a sum over topologies, 
^ complementing the sum over geometries provided by the spin foam model (gravity path integral) alone. Group field 
theories, as argued in 0, |3,0> have a potential and scope that go far beyond that of spin foam models per se, that 
make them fascinating, and that have not been explored but to a very limited extent. 

The coupling of matter fields to quantum gravity can also be described in the spin foam setting 
particular, in the 3d case a precise description of quantum gravity coupled to poiiit p articles of any mass and spin was 
provided in il^i, ll3|, lU , as a spin foam model, based on the classic work of [TsllKjllr^ . and in 0,[3] in the group field 
theory formalism. The basic idea behind the construction of ^2] is to treat Feynman graphs of a generic matter field 
theory as quantum gravity observables, and to obtain a description of the interaction of matter degrees of freedom 
and quantum gravity ones for each line of propagation from the analysis of the classical coupling of gravity and point 
particles, with the latter treated as topological defects of the former. This led not only to a definition of a spin foam 
model (and later a group field theory) for 3d quantum gravity coupled to matter, but also to the derivation of an effec- 
tive non-commutative field theory for scalar fields 14] (and higher spin ones 20]), whose Feynm a.n g raphs amplitudes 
are given by the same quantum amplitudes of the coupled spin foam model, suitably re- written [ij]. The amplitudes 
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for the matter Feynman graphs obtained from the couphng with quantum gravity are, however, a bit pecuhar, in that 
they make use of the Hadamard propagator, instead of the Feynman one, to describe the propagation of matter, i.e. 
they assign the same weights to the external and internal lines of the Feynman graphs, and a rather ad hoc procedure 
is needed to turn these into the more conventional form 14, •211. This can be interpreted, as we are going to prove, as 
a result of a missing causality restriction on the quantum gravity amplitudes one starts with, and that in turn leads 
to a non-causal propagation of matter degrees of freedom at the quantum level, once these have been coupled to gravity. 

Motivated partially by this interpretation of the results of and mainly by general considerations on the role 
of causahty in quantum gravity, and in spin foam models in particular (articulated in H, li^lH,!^), we set out to 
analyze thoroughly the issue of causality in the context of 3d spin foam quantum gravity, with and without matter 
coupling. While the issue we deal with is indeed very general, and we expect the ideas and techniques we use to 
deal with it to be general as well, the 3d setting leads to a particularly neat application of them, as we are going to see. 

Before presenting the details of our analysis and results, we feel it is appropriate to discuss first the role of causality 
in quantum gravity in the path integral formalism and spin foam models, and the general arguments for the need for 
a causality restriction on the histories summed over^, and then the strategy for the coupling of matter fields to 3d 
spin foam models used in |l2lll4| . This will clarify the context, the underlying ideas and the motivations of our work. 



The path integral formulation of quantum gravity entails the construction of a partition function of the form: 



with S{g) being the gravitational Einstein-Hilbert action, for example, for a given spacetime manifold M , and the 
integral being of course over all Lorentzian geometries (metrics up to diffeomorphism equivalence) compatible with 
the given topology. The resulting quantity is in general a complex number and changes into its complex conjugate 
under reversal of the spacetime orientation ^ l27j . It also allows in principle for the computation of quantum gravity 
transition amplitudes (using a canonical gravity terminology) 



where hi and ft,2 are 3d geometries induced by all 4d geometries g on the two (disjoint) boundaries of M (corre- 
sponding to possible canonical quantum gravity states), the complex nature and symmetry properties of the above 
quantities under orientation reversal are the same as those of the partition function. It is obvious that, once defined 
more rigorously, the partition function and transition amplitudes would encode the full dynamics of the classical 
gravity theory at the quantum level, for any choice of spacetime topology M. Also, the above transition amplitudes 
would appear, in a third quantized formulation involving spatial topology changing processes 0, as the 2-point 
function (analogue of the Feynman propagator for a single particle) associated to internal lines of any 3rd quantized 
Feynman graph, confirming that it correctly captures the quantum gravity propagating degrees of freedom (in the 
cosmological context of 9] representing the degrees of freedom of the whole 'quantum universe'). 

A sum-over-histories formulation of the dynamics is possible also in canonical quantum gravity framework. The 
dynamical content of the theory is encoded in the Hamiltonian constraint operator Ti. and one relevant object to 
define for computing physical probability amplitudes, including transition amplitudes, is the inner product between 
physical states, i.e. solutions of the Wheeler-DeWitt equation 7i | = 0. One possible way to define such an 
inner product is in terms of matrix elements of a projector operator from kinematical states onto physical states 
(see e.g. |2^ and references therein), to be thought of as some regularized version of V — J VTe^^^ = S{Ti). The 
sum-over-histories formalism may be used to provide a covariant definition of these matrix elements (including the 



Causality restrictions, one of which is the one on the lapse/proper time/orientation we refer to here, in a discrete path integral for 
quantum gravity, together with the consequences they have for the continuum limit of the same, were shown to be crucial also in 
another approach to discrete quantum gravity, namely Causal Dynamical Triangulations ,25j . 



A. Causality in quantum gravity and spin foam models 
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'vacuum- vacuum'one, equivalent to the partition function), and a covariant definition of physical wave functions 
(satisfying all the quantum gravity constraints, including the Hamiltonian one). The Lagrangian path integral given 
above has then an Hamiltonian counterpart in the case of S x R topology, i.e. when a canonical decomposition of the 
variables is possible. It is given by 



G{h%,hlA = / VNVN 



The (formally defined) integrals are over 3-metrics /ly, their conjugate momenta tt*-' , over the shift vector N''' (this 
integration amounts to the imposition of the spatial diffeomorphism constraint) and over the lapse function N . The 
choice of integration range over the lapse is crucial: the choice of the full infinite range (—00, +00) provides us with the 
result we were aiming for, i.e. a covariant definition of the physical inner product between canonical quantum gravity 
states and thus of the projector operator. In fact the resulting quantity Gjjil^ph]^) = (/12 | hi) ~ (/12 | S{H) \ hi) is 
a solution of the Hamiltonian constraint equation in both its arguments |^ |27| , and it is a real quantity, as expected 
from the canonical inner product. However, Gnih^jjhjj) does not correspond to the lagrangian path integral 
the reason is j23| that positive and negative lapses correspond to the same class of 4-geometries, and the difference 
between the two half-ranges (0, -|-cxi) and (— cxd,0) is only that they correspond to opposite spacetime orientations. 
In other words, the above quantity in Lagrangian formulation is given by a different symmetric choice of quantum 
amplitude, but corresponds to the same set of 4d geometries being integrated over: 



GHihlM^)^ J Vg (e 



-iS{g) 



(2) 



A second choice is possible: one may choose to integrate only in the range (0, -l-oo) to define: 



Gpihjj, hlj ] 



VN / VN' 



Now, the resulting quantity G_f(/ii,/i2) is indeed equal to the Lagrangian path integral and so it also captures 
the full dynamical content of the theory; however, it does not satisfy the Hamiltonian constraint equation in 
either of its arguments; it is instead 0, l27j a Green function for the Hamiltonian constraint operator, satisfying 
7iGi?(^i, /i2) = S{hi — /i2)- It defines thus a quantity analogous to (/12 | j{ \ hi). 

This apparently puzzling situation can be easily understood as arising from the difference between Lagrangian and 
Hamiltonian symmetries, the first being the 4d spacetime diffeomorphisms and the second being the transformations 
generated by the canonical operators Hi and Ti. . The second set of symmetries is actually larger than the first 
and in particular the range N e (0, +00) is perfectly symmetric under transformations of the lapse corresponding to 
4d spacetime diffeos, while it is not under canonical symmetries that can connect positive and negative lapses, thus 
requiring a symmetric range N G (—00, -l-oo). Therefore, while the Lagrangian path integral Z{hi, /12) = Gpih^j, h\^) 
is fully symmetric under spacetime diffeos, it is not symmetric under the transformations generated by the Hamiltonian 
constraint and one needs a further symmetrization Q to satisfy the Wheeler-DeWitt equation. 

Another way to see what is going on is to realize that all this has a very precise analogue in the sum-ovcr-histories 
formulation of the dynamics of a relativistic particle, with G//(/ifj-, /i^^) corresponding to the Hadamard propagator, 
given in momentum space by (5(p^ -|- m^), thus imposing the Hamiltonian constraint equation p^ -\- — 0, and 
Gp^hfj, hjj) corresponding instead to the Feynman propagator, given in momentum space by , thus relaxing 

the same Hamiltonian constraint at the quantum level. 

The restriction to positive lapses is a causality restriction in that it corresponds to imposing that the final 
3-geometry /12 lies in the causal future of the first hi, and imposes a 'timeless ordering' between initial and final 
boundary data; again, it has a precise analogue in the relativistic particle case, where it corresponds to a restriction 
on the particle proper time to positive values only. We want to stress once more another consequence of the restriction, 
which is the consequent relaxation of the Hamiltonian constraint, or in other words of the classical dynamical condition 
in gong over to the quantum regime. 

The Hamiltonian constraint fully encodes the dynamical content of the classical Einstein's equations for the grav- 
itational field (just as it encodes the Klein-Gordon equations of motion for the relativistic particle). The need to 
relax this classical constraint is then understood as the need to allow for purely quantum histories, i.e. for virtual 
propagation. Another way to put it is to say that quantum dynamics takes place off-shell with respect to the classical 
constraint. This is well-known of course in quantum theories of matter dynamics, but it seems not to have been fully 
realized in the context of quantum gravity dynamics. This is true, however, and accepted for 'interacting systems', 
as in quantum field theory. The relevance of these considerations for quantum gravity may be less obvious. There 
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are, however, several reasons why we think 'off-shell' dynamics has to be sought for also in quantum gravity, in other 
words why the Hamiltonian constraint, and the corresponding symmetry, of canonical quantum gravity should not be 
imposed necessarily on quantum gravity states and why the dynamics should be allowed to take place 'off-shell' with 
respect to it. Without entering into much details, let us summarize them. 

One reason is the above mentioned wish to impose causality restrictions, or equivalently a notion of causal, timeless 
ordering among quantum gravity states in computing physical transition amplitudes. Another reason is that, as we 
have seen, such a restriction from the Hamiltonian perspective is necessary to construct a true sum-over-histories 
formulation of quantum gravity dynamics from the Lagrangian point of view, i.e. a Lagrangian path integral, fully 
diffeo- invariant, but with no extra symmetry beingimposed. Also, the inclusion in the formalism of spatial topology 
change by going to a 'third quantized' framework [j, Q involves indeed the inclusion of interactions among universes 
and here is where the covariant off-shell path integral appears as a weight for quantum gravity 'internal lines' of 
propagation, again playing the role the Feynman propagator plays in a usual QFT. One may of course question the 
physical content of such a quantum cosmological framework, but even these doubts are solved by passing to the spin 
foam formalism. In fact, on one hand spin foam models are understood as providing a more rigorous definition of the 
path integral framework for quantum gravity, and thus all the previous considerations regarding causality restrictions 
and off-shell propagation apply to them as well. On the other hand they also accommodate rather easily the presence 
of timelike boundaries, thus providing a local realization of the sum-over-histories idea and thus evade worries about 
the physical interpretation of the quantum cosmology setting. Finally, the group field theory approach and the role 
of spin foam models within it realize also the 'third quantization' idea in a discrete and local way, thus calling once 
more for the need for off-shell propagation of quantum gravity degrees of freedom and for causal spin foam models, by 
clarifying in which sense spin foam models deal with 'interacting gravity degrees of freedom', or better put, interacting 
building blocks of spacetime, that must thus be allowed to propagate off-shell. 

Motivated by all of the above, the construction of such causal spin foam models was initiated in |22j |. continued and 
refined in 23] and then in the group field theory formalism in 24] . The starting point is the realization that current 
spin foam models, including the Ponzano-Regge model for 3d quantum gravity, possess the same symmetry under 
reversal of spacetime orientation that characterizes covariant realizations of the projector operator onto solutions of 
the Hamiltonian constraint. This explains their providing real amplitudes, as opposed to the complex ones one would 
expect from a Lagrangian path integral. All this is fine, and indeed spin foam models were first constructed exactly 
as a way to define in a covariant language the dynamics of canonical loop quantum gravity, but, as we said, there 
are several reasons why one needs to go beyond this initial, albeit important, goal. The mentioned symmetry, as 
appropriate given the discrete setting at the roots of spin foam models, is realized at the level of each dual face of the 
spin foam 2-complex, as a symmetrization between the two opposite orientations it can be given. This was noticed 
in [22] for the 4d Barrett-Crane model together with a consistent way of breaking such symmetry and constructing a 
possible causal version of the same. This breaking was refined to a more elegant and general procedure based on the 
particle analogy, that could be applied to any spin foam model in [23], and finally in a generalized group field 
theory formalism was constructed, that allows for the derivation of causal spin foam models, as well as the usual ones, 
in any dimension and signature, and for a variety of possible quantum spin foam amplitudes. In this process, however, 
the link with the usual path integral formalism for quantum gravity (continuum and discrete) became less manifest. 
In this paper, we derive a causal spin foam model for 3d quantum gravity, that represents the causal counterpart of 
the usual Ponzano-Regge model, starting directly from a discretization of the continuum path integral for gravity in 
first order formalism, and providing a clear discrete counterpart of the continuum causal restriction on the lapse. 

B. Causality in QG coupled to matter: on-shell versus off-shell propagation of matter fields 

One additional motivation for imposing a causality restriction on quantum gravity variables is that such a 
restriction may be in fact needed for describing the quantum off-shell causal propagation of matter degrees of 
freedom of interacting field theories, once these have been coupled to quantum gravity ones and more general spin 
foam models have been constructed. Recall that one example of this is represented by the recently constructed 
coupled Ponzano-Regge spin foam model, as we have mentioned and that we are going to review briefly in the 
following. There are several reason one may expect this at a somehow heuristic level. First of all the dynamics of, 
say, a relativistic particle coupled in the Hamiltonian context is encoded in the coupled Hamiltonian constraint for 
gravity and particle, and it is hard to imagine a procedure (and a motivation for it) that would relax the particle 
part of it, producing the off-shell propagation coded in the Feynman propagator, without affecting and relaxing 
constraints on the the pure gravity sector. Second, the restriction on the particle side that gives the Feynman 
propagator as opposed to the Hadamard one is a restriction on the proper time of the particle, i.e. the spacetime 
length of the particle trajectory, which is intrinsically a gravitational variable/observable, and thus the imposition of 
such a restriction is best interpreted as a restriction on the range of integration of some gravity variables. Even more 



5 



generally speaking, the very notion of time ordering or causal ordering that is at the root of the definition of the 
Feynman propagator for a relativistic particle, and of the time-ordered n-point functions of quantum fields, is in fact 
a quantum gravity notion, so it must find its origin in the appropriate causality restriction on the quantum gravity side. 

Once convinced of these rather heuristic motivations, one naturally comes to expect that, conversely, a missing 
causality restriction on the quantum gravity side would imply an 'a-causal'and ultimately wrong treatment of particle 
dynamics, in particular an a-causal propagation of matter degrees of freedom, once these have been coupled to 
quantum gravity in some spin foam model. From this perspective, then, an apparently puzzling and a rather disturbing 
feature of the coupled Ponzano-Regge model of [Tj, 13 looks instead natural and confirms the above motivations and 
expectations. The feature we refer to is the fact that, as we are going to see in a few more details in the next 
section, the Feynman graphs for the effective field theory describing the matter degrees of freedom have Hadamard 
propagators instead of Feynman ones weighting their internal lines, and thus describe these degrees of freedom as 
being always on-shell, i.e. always solving the classical Hamiltonian constraint/equations of motion. This also implies 
a sort a strong a-causality in the sense that no ordering of their arguments is encoded in the Feynman amplitudes 
(and correspondingly there is no difference between particles and anti-particles). This turns out to be true both at the 
non-commutative field theory level, encoding the full quantum gravity corrections to matter propagation, and in the 
no-gravity limit, when one indeed recovers a matter field theory in flat space but, as said, with Hadamard propagators 
in place of Feynman ones. This means that the result cannot be attributed to any sort of quantum gravity corrections 
to ordinary matter propagation, but it has really its origin in some more basic feature of the formalism. 

Of course, having the Hadamard propagator on each internal line, one can try to reconstruct from it the corre- 
sponding Feynman propagator, and this is indeed done in [T3.l2l||. However, on the one hand this is necessarily a very 
ad hoc procedure and does not teach us much about the origin of this puzzling feature of the model nor about its 
meaning, and on the other hand it is not a unique procedure in that many different functions representing the non- 
commutative version of the usual Feynman propagator can have the same on-shell (Hadamard) counterpart, because 
a generic p-dependent factor may become a constant function of the mass m only, due to the — constraint 
being imposed, as we are going to see and discuss, and because such p-dependent factors may also disappear in the 
no-gravity limit, and thus such different non-commutative versions of the Feynman propagator may end up having 
the same flat space counterpart. 

We will see in detail how our ah initio procedure solves all these ambiguities. The usual coupling of particles to 
3d quantum gravity, plus the same causality restriction of the pure gravity case, automatically provides us with a 
well-defined causal version of the coupled Ponzano-Regge model, and with causal Feynman amplitudes for matter 
fields, with a unique prescription for the non-commutative version of the Feynman propagator, and one that of course 
has the right and expected no-gravity limit. 



Having presented the motivation for our work, let us now outline our results as we are going to present them. After 
a very brief review of the spin foam formulation of 3d quantum gravity coupled to matter fields (section II) , we will: 
1) identify the correct procedure for implementing a causality restriction in the spin foam formulation of pure 3d 
quantum gravity (section III); 2) construct explicitly the resulting causa/ spin foam model (section HI); 3) add matter 
degrees of freedom to it and construct a causal spin foam model for 3d quantum gravity coupled to matter fields^ 
(section IV); 4) show that the corresponding spin foam amplitudes admit a natural approximation as the Feynman 
amplitudes of a non-commutative field theory, very similar to that obtained in |l4j . with the appropriate Feynman 
propagators being used (section V); 5) show that this effective field theory, and its Feynman amplitudes, reduce to 
the usual QFT action and Feynman amplitudes in flat space, as they should, in the "no gravity limit" (section V). We 
then conclude with some outlook on possible developments of the results we have presented. 

II. SPIN FOAM FORMULATION OF 3D QUANTUM GRAVITY COUPLED TO MATTER FIELDS 

Classical 3-dimensional gravity can be formulated in first order formalism as a BF theory. 



^ We deal explicitly with scalar matter only, but none of our procedures or results is significantly affected by the presence of spin degrees of 
freedom, which simply modifies some of the details. Also, we work in the Riemannian quantum gravity setting, exclusively for simplicity 
of notation and of mathematical details. We stress however that none of our results or procedures used is affected by a switch to the 
Lorentzian signature, that causes only a slightly greater complexity of the mathematical details, ultimately due to the non-compactness 
of the relevant group, i.e. SL{2, R). 
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The dynamical variables of the theory are: 

• an SU (2) connection, denoted by A, which can be locally thought as an su (2)-valued 1-form; 

• an su(2)-valued 1-form, denoted by E^, 
and the action is given by: 

S{E,A,M)^ [ Tr{E/\F{A)) (3) 

for a generic spacetime manifold M, where F{A) is the curvature of A, thus a Lie algebra-valued 2-form, and the trace 
is taken in the adjoint representation. The theory is quantized using the path integral method, i.e. by the partition 
function: 



The spin foam representation of the theory is obtained by a simple discretization of the above classical and quantum 
theory 0, of course, the discrete setting helps to a great extent in making mathematical sense of the above 

functional integrals. So, we introduce a triangulation A, and its topological dual A*. Keeping in mind that the 
geometric information of a simplicial complex is fully captured by the set of edge lengths, that the curvature is 
located on (d-2)-simplices (thus again along edges in 3d), and that £' is a 1-form, while F is a 2-form, a natural way 
to discretize the variables is to do the following: 

• integrate the E field along the edges of the triangulation, obtaining su (2) elements, one for each edge. 

• integrate the F field over the dual faces (i.e. over the faces of the dual cellular complex), obtaining an su(2) 
element for every dual face (dual to an edge of the triangulation) . 

Given that, under any lattice refining, (^Holonomy = 1 + Curvature + . . . ^ , a convenient but also appropriate choice 

is to define the discretized curvature F through the holonomy around the boundary of the same dual face on which 
it is discretized, which in turn can be obtained from a product of parallel transports (path-ordered exponentials of 
integrals) of the fundamental connection A along dual edges, that give in the end SU{2) group elements. 
Thus in the end a natural set of variables for the discretized theory is: 

• an su(2) element for every edge e, denoted by X^^. 

• an SU(2) element for every dual edge e*, denoted by t/g. . The product of all the g^* 's corresponding to the same 
e in cyclic order, i.e. the holonomy around a dual face dual to the edge e, will be denoted by Gg. 

With this discretization the action above reads 

5=^Tr(XeGe) (4) 

edges 

Now, quantization is straightforward. The partition function is given by^ 

Z = n /rf^e n / d9e> (5) 

The integral over in lO can be performed^, to give: 

Z = lll dg,.Y[5{Ge) (6) 



^ This is true only locally. Globally it is a section of the bundle associated to the principal SU(2) bundle (over which A is defined), via 
the adjoint representation. 

The factor is needed in order to get a delta function over SU(2). This point was emphasized also in 12j, to which we refer for 

more details. e(G) = signiC os{9)) where 9 parametrizes G (see Appendix A). 
^ Note that this is true so easily only in the pure partition function, i.e. not in general if one is computing a generic observable that is a 
function of the as well as of the g^* . 
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Now the delta functions on the group can be expanded in irreducible representations using harmonic analysis on 
SU(2), and then the integrals over the g^' can be performed as well, leaving us with the spin foam version of the 
same partition function |l2| : 



z=y IT ( 1 



{jg} e tctrahcdra 

where the quantity 2je + 1 is the dimension of the representation je labeling the edge e of the triangulation, and 
the weight for each tetrahedron in the triangulation is given by the 6_;/-symbol, a function of the six representations 
labeling the six edges of the tetrahedron itself. This is the so-called Ponzano-Regge spin foam model. 

The coupling of matter fields is obtained, as anticipated, by treating a full Feynman graph F of a scalar field theory, 
with its hidden dependence on geometric variables, as a quantum gravity observable. The coupling between geometric 
and matter degrees of freedom at each line of propagation of the graph is obtained by a discretization of the continuum 
action describing the coupling of gravity to relativistic point particles in 3d, with the line of the Feynman graph thus 
being interpreted as the trajectory of a relativistic particle in a 3d spacetime, and by the subsequent integration over 
particle data. 

One therefore introduces curves 7i's embedded in the spacetime manifold M, where 7i is interpreted to be the 
worldlinc of the i-th particle. Collectively, these curves make a 1-dimensional graph F, interpreted as a Feynman 
graph of a field theory. F will often be called the 'particle graph' below. The analogue of JSJ is then^ 

S^l Tr{EAF)+y2rni I Tr[{E + dAq) uJqu-^] (7) 

where, q is an su (2)-valued 0-form (function) and u is an SU(2)-valued 0-form. miuJou~^ is interpreted as the 
classical (thus already on-shcU) momentum of the i-th particle and thus shall be denoted by pi'^ . The equations of 
motion, which are obtained after varying the E and the q are 

FiAix)) = ^p,<5(a;-a:^J (8) 

i 

dAPr = (9) 

The first equation means that the curvature is flat everywhere except at the locations of the particles where it has 
a conical singularity. The second equation means that the momentum of a particle is covariantly constant along its 
worldline. 

The full continuum partition function is now 



j VAVEY[{Vq,Vp,) 



Note that in the path integral we have something of the form 



The effect of the integral over g^'s is then simply to impose equations However, as will be apparent in a moment, 
the discretization itself, if chosen appropriately, will impose these constraints.^ Thus we shall drop the last term in 
Q and discretize a slightly different but equivalent theory given by: 

S = [ Tr{E/\F)~y] [ Tr{Ep{) (10) 
Jm ■ Jii 



^ This action was first derived in llTl . for a nice discussion of it see the beginning of '29'|. 

rhi = AizGrrti, where rui is the mass of the i-th particle and G is the Newton's constant. We shall use in this section units in which 

4ttG = 1, and hence will call rhi the mass as well, while it is in fact the deficit angle associated to the mass rUi. We shall restore 

Newton's constant to the equations in section V since we will be interested in the limit G — > 0. 
* Essentially, this is because each edge of the triangulation (adapted to the curves, or more precisely to the Feynman graph) will carry a 

single group element, to be considered the discretized analogue of p. 
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To discretize this action we, as before, introduce a triangulation A and its dual A*. We insist, however, on 
a triangulation that is adapted to F, i.e. the 7i's coincide with some of the edges of A. The discretizations of 
the E and A fields are given as before. The discrete analogue of pi is given by u 6^™''^°^ u^^, the group element 
corresponding to the su (2) element [fhiU Jqu~^) It is important to note that there is one u for each edge 

of r, i.e. the discretized pi is constant along each edge of F, thus satisfying l^. Since there is exactly one 
associated to every edge of F, from here on we shall label the mass of the particle by the edge that it is propagating on. 

With these definitions, the analogue of Q for action (fTUI) is given by^ 

S= Yl Tr{X,Ge)+ Tr{X,GeU-^e'^^-^°^uJ (11) 

edges edges 

It is obvious that if one calculates the equations of motion associated to the Xe variables on the particle graph one 
obtains Ge = u^e^'^^'-'°^u~^ = Uehi,u~^, which is essentially the holonomy version of (jSJ. 

The partition function for the coupled theory is analogous to (jSJ: 

e e* e ^ eeF 

Also in this case the integral over Xg can be performed explicitly (again, for observables that are not functions of 
them) to give: 

Z = l[ f dge' n / n C^-) n ^ (GeUekeU-^) , (13) 

e' eer-' e^F e6F 

and also in this case one can obtain a spin foam presentation of the same partition function using the harmonic 
analysis on SU(2) and thus expanding the delta functions in H13() in irreps: 



{je} e^T eGF tctrahod 

This is the so-called coupled Ponzano-Regge model the case of massive spinless point particles, modulo some 

arbitrary (at this stage) normalization factor given by a function of the mass variables TOg {l2l |Q| . 

One key result, obtained in (lj|, is that, having re- interpreted the G^ at the location of the particles as their gravity- 
induced physical momenta at the quantum level, the above partition function can be re-expressed as the Feynman 
evaluation of a Feynman graph for a non-commutative field theory. The way this is achieved is by introducing non- 
commutative momenta valued in Bi^i^^), i.e. the 3-ball in M^" with radius 1/k = 47rG, through: g = y'l — K^|_Ppl + 

kP ■ J, the star product ea^^''^^^^^-' ★ e^'^^'^^^^^ — e^'^^'-^^^^^^ and the 'Fourier' transform naturally associated to 
the above product f{X) = J dg e~^'^^'^-^^'^^^^ f{P{g)). The result (modulo pre-factors) is: 

/ ^^HMI^^I ^ ^^{®e<,O.Pe). 

eGF"^ eGF ^ ' uGF 

Basically, the delta function enforcing the conservation of the non-commutative momenta at the vertices of the 
Feynman graph results from the delta functions for pure gravity edges of the triangulation in p3|l , while the Hadamard 
propagators for the lines of the particle Feynman graph result from the delta function coupling gravity and matter 
degrees of freedom in H13|l . The problem we had anticipated is now manifest. Hadamard propagators are not the 
correct way to capture the quantum dynamics of matter degrees of freedom, in particular they possess the wrong 
causal properties, and are not what comes out of the perturbative expansion of a QFT for matter fields. The problem 
is, as we said, bypassed in |14L I2H by re-constructing a possible definition of the Feynman propagator from the 



The addition of su(2) elements F and Pi is converted into product of SU(2) elements Ge and Ug e 



(me Jo), 
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obtained Hadamard propagator, and using this possible definition to reconstruct the non-commutative QFT action 
that would originate the above Feynman amplitudes. However, as mentioned in the introduction, the reconstruction 
is highly non-unique, in particular because the definition of the Hadamard propagator is non-unique: the insertion of 
any function of is turned by the delta function into an irrelevant insertion of a function of mg. Also, the choice 
made in [T^ l2l| is natural and gives the correct result in the no-gravity limit, i.e. G ^ 0, but, again, it is clear 
that this limit does not constrain the initial choice enough to achieve uniqueness: any insertion of a function of 
that gives a constant in the limit docs not affect the result. We will see how a correct ah initio derivation of a causal 
coupled model will solve naturally and uniquely all these difficulties. 



III. PURE GRAVITY - CAUSAL PONZANO-REGGE MODEL 



We want now to derive a causal version of the Ponzano-Regge model for pure 3d quantum gravity, starting from 
the discrete formulation of 3d BP theory reviewed above, and imposing on it the causality restriction that reduces 
the corresponding path integral to that of 3d quantum gravity, in Lagrangian form. By deriving everything directly 
from the discretization of the continuum theory, we will avoid any sort of ambiguity (apart those intrinsic to any 
discretization procedure) and obtain a unique prescription for how the causal Ponzano-Regge model should differ 
from the usual one. 

Let us now go back to the partition function: 

Z{M) ^ j VE j PAe^^(^'^'^^), 

and its discrete analogue Q: 



e e* e 

Another fruitful way of looking at it is to consider the action in ||SJ) a discrete analogue of 

J dt J d^x{N'H^ + NH) (15) 

Where N\ N are the shift and the lapse; and Hf, are the diffeo and the Hamiltonian constraints. This is almost 
immediate from 0, since after one integrates over the Xe's one is left with (5(Ge)'s which are essentially imposing 
the constraints — 0, or Hf = = H"^ . This line of reasoning has been explored in detail in ^5], and indeed it 
was verified that the above partition function imposes on the spin network states the diffeomorphism as well as the 
Hamiltonian constraints. 



However, for reasons mentioned above, we do not want to impose the Hamiltonian constraint. Or, put differently, 
we want to obtain the causal analogue of the above partition function. At the Lagrangian level this means restricting 
the geometric variables in the partition function to those corresponding to a positive orientation of spacetime only, 
i.e. constructing a discrete analogue of the constrained partition function: 

Z{M) = JvE J VAQ{det{E)) e'SiE.AM)^ 

where the Q{det{E)) is the Heaviside step function effectively imposing the constraint det{E) > 0, i.e. the positive 
orientation of spacetime as encoded in the triad. Before we perform such a causal modification of the Ponzano-Regge 
model, let us briefly review its symmetries |32j . In the continuum, PR has two symmetries: 

• The usual SU(2) gauge transformations, given locally hy E gEg^^ and A — s- gAg^^ + gdg~^, where g is an 
SU(2) group element. This symmetry will be called 'rotational' symmetry below. 

• An additional symmetry, which can be seen as originating at the classical level from the Bianchi identity, 
E ^ E + dAfp, where </> is an su (2)-valued function. This symmetry will be called a 'translational' symmetry 
below. 
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The above two symmetries combined, are equivalent to difFeomorphism symmetry on shell . 

Let us write X^. = Xe ■ J and Ge ~ Cos{9e)l + fie ■ JSin{9e), where the J's are the su (2) generators and de and He 
are the standard parameters for the SU(2) group. With these definitions, action becomes 

S = — Xe ■ fie Sin{9e) (16) 

edges 

Now, the symmetries become 

• Simultaneous rotations of Xe's and Tie's, i.e. Xe — > OeXe and rie Oefte, where Oe is an orthogonal transforma- 
tion associated to every edge. This is the reason for calling this symmetry rotational. 

• The translational symmetry has a more complicated expression [32j . but essentially it IS X fp ^ c I 0g ■ 

It should be clear from the discussion in the introduction that the way to impose causality is to put some restriction 
on the Xe integrations. This is because we would like to restrict the integration to 'positive' triads, and triads 
correspond to Xe's after discretization. Obviously, a restriction to the positive triads will not affect the rotational 
invariance, since it is impossible to convert a positive triad to a negative one by a rotational gauge transformation. 
Thus, the causal theory has to still be invariant under the rotational symmetry. On the other hand, the translational 
symmetry has to be broken, at least partially, since it is clear that one can go from 'positive' triads to negative ones 
by doing a translational transformation (just choose d-Aff' — — 2_B). 

Another way of looking at it is to say that we want to restrict to positive lapses only. Now, comparing H15|l to 
(|16|l one could identify {N,Ni) with x and {H,Hi) with nSin{9). Thus, the translational symmetry implies sym- 
metry under N N+constant. So, it is evident that there's no way to reconcile this with the fact that we want TV > 0. 

So, how should we restrict the range of Xe's, to implement the restriction detE > 0, while keeping the rotational 
invariance? Before the prescription is provided, note the following four points: 

a - In the continuum, if ii^ — > —E then det{E) —det{E). Hence, if Xe can be thought of as the discretization of a 
positive triad, then —Xe should be considered the discrete analogue of a negative one. Thus, if a certain subset 
of corresponds to positive triads, its parity-inverted image through the origin corresponds to the negative 
ones. 

b - In the continuum, if det{E) > then det{XE) > as well, for any A > 0. Hence, if Xe corresponds to positive 
triads. Ax does so as well. 

c - The rotational invariance implies that if Xe is a positive triad, then Oe x is also a positive one, where Oe is 
any rotation with fie as its axis. Note that [b] and [c] taken together imply that, the region of all x's which 
correspond to positive triads is a union of cylindrical cones with axis n. 

d - Finally consider the continuum action, with a positively oriented triad. Then start 'deforming' the triad until it 
becomes negative. Since, the Lagrangian is proportional to the volume form, this means that the Lagrangian, 
switches sign as E is deformed into —E. This implies that the Lagrangian passes through zero, along the 'path 
of deformation'. This implies that if we start from the discrete Lagrangian and start 'deforming' the x, then as 
long as the lagrangian doesn't vanish, we haven't switched from positive triads to negative ones. The earliest 
we can do so is the moment the Lagrangian passes through zero. 

Thus, assume that x is almost along the corresponding n, and that this x corresponds to a positive triad. Then, 
[b] and [c] imply that all a?'s within a very 'thin' cone with axis n are positive as well. Keeping in mind [d], start 
'opening' the cone. Then, as is evident from H16|) . the earliest the Lagrangian can become zero is when x is orthogonal 
to n. However, recalling [a], one sees that this is the largest allowed region of positive triads (otherwise the regions 
where the non degenerate triads are positive and where they are negative overlap). 

The restriction to positively oriented triads is obtained therefore by considering only those x that make a positive dot 
product with n. In other words, one has to impose that the projection of the triad along the normal to the surface to 



Xe is an element of R'', He is an element of 5^, and 6e is a real parameter which goes between and tt. Check Appendix A for more 
details. 
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which the holonomy defining the curvature refers is positive. Notice that this is exactly analogous (although a proper 
Hamiltonian treatment of the discrete action we are dealing with would take us too far) to the continuum restriction 
on the lapse function, defined indeed as the projection of the triad field along the normal to the surface to which 
canonical variables refer. 

We arrive at the following prescription: 

In the partition function J^J) one has to integrate only over those 's such that Xf. • rig > 



If the above geometric arguments were not sufficient to derive and motivate our prescription for the causal- 
ity/orientation dependence restriction, this could have been motivated also by comparing the action 1)16(1 to the 
Regge action in first order form, in terms of group-theoretic variables, derived in Isij. The expression is: 



S= Lf,Sin{Qe 



(17) 



edges 



Where, Le is the length of the edge e, and 6e is the deficit angle at this edge. Comparing l[T7|l to ((T?)|l . and identifying 
Lf. with Xe • ??e, it is immediate that the above prescription corresponds to integrating only over positively oriented 
edge lengths, exactly as one would expect^^. 

We shall perform the Xe integrals in © . Since the partition function factorizes per edge (dual face) , and gives the 
same weight to each, it is sufficient to do only one of these integrals, for an arbitrary choice of edge, henceforth we 
shall drop the subscript e from our vectors. 

So, let us now compute the integral over x: 



cPx g~«(^-"'S'™(f)) _ 



x-n>0 



g-ix-n(Sm(e)+ie) 
x-n>0 

^-ircos{a){s^nie)+ic) ^2 sin{a)drdadp 
~2m 



Sin{9) + ie 

-2711 

Sin{9) + ie 
-27ri 



( \ L 

\ii{Sin{e)+ie)) \i(Sin(9) 



2 / 1 s2 

ij{Sin{e)+ie)J \i{Sin{e) + ie) 
1 



1 



{Sin{e) + ie)^ \ 72 
Cs I 



4 [Sin{e)+ief 



(18) 



(19) 



In the above 7 is the cosine of the upper limit of the a integration, 7 = Cos( J — 5) — Sin{5), where (5 is a small 
positive number, the factor of a quarter is included to recover the usual normalization used in the Ponzano-Regge 
model if we take the sum of H19|l and its complex conjugate. Due to the importance of the above result, we give a 
more rigorous derivation of it in Appendix B. The regulator e plays exactly the same role as the regulator in the 
Feynman propagator in the usual quantum field theory, as will be apparent in section V. It is understood that it has 
to be sent to zero at the end of every calculation, i.e. when computing normalized expectation values of observables, 
including field correlations. Note that the presence of a small explicit regulator e is fully equivalent to having 9 



I31I considered the SO(3) theory, while we are considering the SU(2) one. This is the reason that 9^ in 1161 is not ©e in I17i . In fact 
©e = 2 6e. Classically the two theories are, of course, the same, since the equations of motion, in both cases, imply that ©e = = 6e. 
Note that g [0, tt], consequently the deficit angle £ [0,27r]. This should be clear, since after parallel transporting a vector around 
the edge, there's no difference between a rotation by angle 6 or by 9 + 2-nn. The fact that in the usual Regge calculus the range of the 
deficit angle is [—00, 2it\ is an outcome of the traditional second order formulation. 
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analytically continued to the complex domain and having acquired a small positive imaginary part. 

In the above S has to be set to 0, but it is clear that this gives an infinity. Fortunately, all the infinity is contained 
in a constant which factors out and does not contribute to the physical calculations wc will be interested in. Let us 
clarify: we will get this constant (which we have denoted by Cs) for every edge of our triangulation, thus the partition 
function (which becomes a function of 5) will be given by 



Zs = (Cs) 



E\ 



where the above equation defines Z, and ji?! is the number of edges in the triangulation. 
Now as far as observables of the connection (A in the continuum, g's in the discrete) are concerned we can 
as well use Z for calculations, which is independent of S. This is because as far as the connection is con- 
cerned the difference between the different i5's is a connection-independent constant and thus will cancel from any 
calculation^^. Since in this paper we will be interested only in connection-related questions, we will use Z from now on. 

What about the observables which are functions of the x's? In that case we of course cannot use the above 
argument and will have to integrate all the way to S — for those x in the argument of the observable. Whether 
the expectation value of such observable converges or not is a question which depends on the observable under 
consideration, and cannot be answered a priori^"^. 

Keeping equation H19|) in mind, the partition function Q becomes 



Z 



11 d9e^ 11 3 11 IJsMeJT^ ^20) 



e*eA* e eSA 



One could expand now the above expression in characters, then proceed in the usual way, by integrating over gg. 's 
and obtaining 6-j symbols. The only difference from the usual Ponzano-Regge partition function now would be in the 
weight assigned to each edge (dual face). So, if the dual face is labeled by a representation j, instead of the dimension 
of the representation one gets more complicated functions of j's and e's, denoted by Cj{e). Consequently, the spinfoam 
version of is given by 

^= E n^.(^) n {6j} 

labellings eEA tetrahedra^A 

The sum goes over all labellings of the faces of the dual two-complex by representations of SU(2). In the above 
{6j} denotes the usual 6-j symbol, and the Cj(e)'s are given by the following formulae 



(l + e(G)) l I 

2 4 {Sin{9) + ie)^ ^ 



3 



Sin{9) ' TT Jo A {Sin{0) + ie)^ 

The coefficients Cj{e) can be given as linear combinations of hypergeometric functions. 
Going back to 12()|l we note the following properties of this expression: 

I The amplitude is sharply peaked around the classical, on-shell, configuration, 6^ = 0. Exactly as one would 
expect. Heuristically, one can think of the edge amplitudes as 'smoothed' or 'regularized' Dirac delta functions, 
just as the Feynman propagator for a spinless particle p'i^„^i can be considered as a smoothed S{p^ — m?)}'^ 



As is well known, for any observable 0(A), the expectation value is given by (0(A)) = •- — ^ pAPEe'^^ — ' '^'^"^ overall connection- 
independent constant will be irrelevant for such calculations. 
The observable should go to zero faster than j-kr^ as Ixl — > oo. 

The Feynman propagator is equal to a delta function + a principal value distribution, which can be thought of as a quantum correction 
to the 'classical' delta. 
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II The amplitude generically is a complex number, moreover Z ^ Z when the orientations of the Xe's are reversed. 
This is is not obvious from H2U|) since the Xg's have been integrated out. However it is clear from the calculation 
on the previous page, that if one replaces x with —x one gets the complex conjugate of the final answer. Again, 
this is exactly what one would expect from the considerations about the continuum path integral, and matches 
the result we were aiming to obtain, as discussed in the introduction. 

III The amplitude is not triangulation independent. Although this may seem surprising, in fact it should be 
expected, since we are not on shell anymore. 

One expects triangulation independence because of the fact that classically there arc no local degrees of freedom, 
however to show that this is indeed the case one needs to invoke the equations of motion (constraints) in order 
to show that there are no local degrees of freedom left after gauging away the symmetries. Since we are not 
imposing any constraints we have local 'virtual' degrees of freedom, even though classically the theory is totally 
trivial. This is just like the case of a relativistic particle. Classically the magnitude of the momentum is fixed 
to be the mass, but quantum mechanically it can be anything. 

Another way of looking at this, and to understand how no modification (and no restriction) of the BF partition 
function, and in particular no realization of the Lagrangian path integral for pure gravity, could be triangulation 
independent, comes from the canonical picture, along the lines introduced in [3. If one tries to write an 
'evolution' operator U, as opposed to a projection operator onto solutions of the Hamiltonian constraint, using 
the same approach^^, then the rotational gauge invariance implies that U has to be expandable in characters. 
However the 'fusion' identity will only be satisfied, as proven in P|, if Cj = 2j + 1, i.e. for the delta function 
over the group, imposing the Hamiltonian constraint, as in the usual Ponzano-Regge model, or Cj = 0, thus 
giving a trivial operator. Without the fusion identity, there can be no regularization independence and hence 
no triangulation independence. 

IV In the above, we have totally glossed over the issue of gauge fixing We have mentioned above that the 
translational symmetry is broken by our restriction on the x variables. However, this is only partially true. 
Roughly, the translations along the n are broken, while the ones in the directions orthogonal to n remain intact. 
Thus, what one should do, is to choose a maximal tree T, choose a gauge fixing condition to be implemented 
along this tree and compute that Faddeev-Popov (FP) determinant. If we choose the simplest gauge condition, 
which is {x — X ■ n — Q\ i.e. we set the the component of x orthogonal to n to zero, then using the results of 
[3^ it is easy to show that the FP determinant has the form A = neGT(l + l^eP), where fig is a complicated 
function of all the Gg's associated to the edge e^^. In this choice of gauge fixing, the amplitudes for gauge fixed 
edges are given by a regularized g^^^ instead of the cubic inverse power of the un-gauge-fixed amplitudes, as 
it is easy to verify by direct calculation. 

It can be shown that in the case when we have delta functions everywhere as edge amplitudes, i.e. by using the 
classical equation of motion, then the FP determinant evaluates to one. In our case this is no longer the case, 
since our edge amplitudes have support away from the flat connections. Since the FP determinant is a nonlocal 
function (what we mean by this is that it cannot be factorized into a product of contributions each of which 
is a function of variables coming from a single edge), it depends on the triangulation and does not simplify. 
The calculation of this Faddeev-Popov determinant, and thus the presentation of the fully gauge fixed causal 
partition function, will be left for future work. 



IV. MATTER COUPLING - CAUSAL COUPLED PONZANO-REGGE MODEL 

We now consider the case of matter coupling, and show that doing a restriction similar to the one done above 
one obtains a model similar to the one proposed in |l4| . but with the correct Feynman propagators instead of the 
Hadamard ones obtained there. 

We shall now couple gravity to matter along the lines introduced in |12| . and reviewed in section II. So we end up 
with the partition function: 



Essentially, one must define some regularized version of {\U\) = (| / dre""^!) = (| jj-^^^ I); where H is the Hamiltonian constraint. 
^® Effectively it is given by Baker-Campbell-Hausdorff formula. What we mea n by 'the Ge's associated' is the set of all Ge's which 
co rres pond to all the edges meeting at the vertex of the original edge. See l32l for details. The only difference from the formulas given 
in I32I is that our f2e is the projection of the one defined there on the plane orthogonal to associated to the edge e. 
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Once more, before showing what type of causahty constraint we are going to impose on this theory and the result 
of this imposition, let us discuss very briefly the symmetries of the above theory, because it will help in clarifying 
some aspects of the causal restriction and of the resulting causal model. 

The theory retains the symmetries of the pure gravity case, i.e. the rotational and translational ones, if one augments 
the previous transformations with: 

• (/—!■ gqg^^, and u — )■ ug^^ , or equivalently p —> gpg^^, for the rotational symmetry. 

• q ^ q — 4>, for the translational one. 

Before we give the partition function, let us write our variables, analogously to the previous section, as: 

Xe = Xe-J (21) 

Ge = Cos{e)l + He ■ JSin{e) (22) 
^-ig(A,Jo)y^ = Cos{me)l+Ue- JSiniihe) (23) 

With these definitions, the action Hll|) becomes 

S = — ^^-Ee • neSin{9e) — ^Cos{rhe)Sin{0e)Xe ■ + Sin{rhe)Cos{Oe)Xe ■ Ue — Sin{ifie)Sin{6f.)xe ■ {fie X Ue) 



Spurs gravity — Sp S gravity coupled to matter — Sc 

(24) 



The partition function is 



n / rf.9e-n^^^y^n / ^-^^n / ^ 



Where, 



Zp ^ \ \ I dx,, e*^" 



e^r-^ 



The above separation of the partition function is possible because the action (I24f) is factorized per edge (dual face) 
as far as the Xe and Ue variables are concerned. We want to do the integrals over Xe's and the -Ue's. The treatment 
of Zp is exactly the same as it was before, and was shown in the previous section. 

We move on to Zc- Now Zc can be written as a product of contributions coming from each edge of F. 

^e=n^e 

eer 

where 

Z":. = 



where is the part of the action H24|l coming from edge e. We shall concentrate on only one edge of F and, 
consequently, drop the 'e' labels on the variables from the formulas below. 
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^ — i^Cos(fh)Sin(6)x-n+Sin{7h)Cos(9)x-u — Sin{m (25) 

We want to do the integrals over x and u now. What are the ranges of integration? For u the answer is simple: it 
is S"^ ~ SU{2)/U{1), where the C/(l) chosen is the subgroup of SU{2) to which h = e™-^" belongs. What should the 
range of integration for x be? From considerations in the previous section, we expect it to be half of M'^, the question 
is, which half? 

We know that, in the continuum, the component of E along the trajectory of the particle, is equal to the proper 
time. This means that the component of E along the momentum of the particle gives the proper time, and it is 
this component that we must restrict to positive values. We have two potential vectors to call 'momentum', F and 
p (see (llOll '). At the classical level, there is no difference between F and p, since the equations of motion set 
them equal to each other. It might seem at first glance that p is the right momentum to use. However, it turns 
out that, since we are interested in the quantum theory, F is the right choice. This is because the magnitude of 
p is fixed to be m, i.e. p is already on shell. But we know that in quantum mechanics the momentum of the 
particle can be anything it wants and is not restricted to the classical value. Thus, we conclude that F is the more 
suitable candidate for momentum. Of course, this is confirmed by the whole body of knowledge about quantum par- 
ticles coupled to 3d quantum gravity, as can be seen in [l2. Ila lla Il7l | . as well as by the results we are about to present. 

It follows from the above that we must restrict the component of E along the F to positive values only. Going 
through the same arguments as in the previous section, it follows that in the discrete case this means restricting the 
integration in the partition function in such a way so that x ■ n > 0. Notice that this is exactly the same prescription 
that we have used in the pure quantum gravity case. This confirms once more that there is a universal (i.e. valid 
in presence as well as in absence of matter) restriction on the geometric variables that produces the correct causal 
quantum amplitudes we are after. 

The upshot of the above is that we do the x integral over the half space for which x ■ ft > 0, remembering^'' that 
this implies that the n has to be complexified, or equivalently that the z — n ■ u is complex (hence the suggestive 
notation). The result of integrating H25|l over this range of x gives^^ 

where 

62 = ^(ri ■ ufSin^{9)Sin^{rh) + 2{n ■ u) Sin{6)Cos{e)Sin(m)Cos{rh) + 1 - Cos^{m)Cos^{e). 
Substituting this expression for b into (|26|l . and using spherical coordinates to do one angular integral we get 



l+ie 

dz i 

-l+ic 



y/~Sin^{9)Sin^{m)z^ + 2Sin{d)Cos{9)Sin{m)Cos{m)z + 1 - Cos^{m)Cos^{9) (27) 
Before we analyze this integral, let us simplify it somewhat. So, the integral above is equal to 



-l+ie 



1 - Siv?{e)Si'n?{m) z - Cot{e)Cot{m) 



2\ 2 



-Cos(9+m)+it 

'iSin{e)Sin{fh) J -cos{e-fn)+te 



dw 



1 — w'' 



(28) 



Once more, it is crucial to keep in mind that the above is a complex integral, otherwise the step from (|25|) to (|26|) is 
simply incorrect. This comes from the need to complexify the group element (its parameters, and thus n) in order to 



See Appendix B. 

See Appendix B again. We have implicitly divided by the constant Cs from the previous section. Again, the factor of half is needed if 
one wants to recover the coupled Ponzano-Regge with usual normalization for the amplitudes (no numerical constants in front of the 
delta function). 
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make sense of the integral over x that gives rise to (|27() . In other words, we are not trying to evaluate a real integral 
using complex methods. The integral is complex to start with, and forcing it into the real domain creates numerous 
problems, which are simply not there if one deals with it correctly. The main problem with considering the integral as 
it were in the real domain is that one loses analyticity: the real function resulting from performing (|27l) as it was a real 
integral is not only discontinuous, which is normal as it has poles, but it does not admit any analytical continuation 
even away from the poles. As such, it cannot be considered as a propagator (Feynman or else) of any field theory. 
In fact, as a consequence of the non-analyticity, the action of the theory, in momentum space, which would have it 
as its free propagator, would be discontinuous across the point where the momentum is equal to the mass (classical 
solution of the field equations) . More precisely, the action switches from the one for a usual scalar field theory when 
the momentum is less than the mass, to an action of a nonlocal theory when the momentum is larger than the mass. 
All these problems persist even in the no-gravity (G —^ 0) limit. 

All these issues are solved automatically once one realizes that (|26|l is a complex integral from the very beginning. 
Once this has been done, one has then to specify the contour of integration, that will depend on the two parameters 
9 and to. Moreover, the ambiguity in our case is even greater, since the integrand is a complex function with branch 
points, and the natural domain of integration becomes the associated Riemann surface and not the complex plane. 
Thus there are even more possible paths of integration in this case, since we can go from one sheet to the other. 

The upshot of the above is that after taking analyticity into account there are essentially only two choices of the 
contours of integration. One choice gives the usual Feynman propagator as we shall see. The other seems to give 
the difference of advanced and retarded Green functions (instead of the sum, which is the Feynman propagator), and 
would correspond to a non-local field theory even in the abelian limit. We show the contours of integration on the 
Riemann surface and on the complex plane in Appendix C. 

If one follows the contours prescribed, one obtains as a result of the z integral, 

Cos{e) i 



2Sin{m) Sin{9 — to) Sin{6 + m) + ie 
Or, after trigonometric manipulations 

Cos{9) i 



(29) 



2Sin{m) Siii?{9) — Si'n?{rh) -f ie 

This last expression, agrees, up to a Cos{9) factor which is constant in the abelian limit, with the expression for 
the nonabelian Feynman propagator put forward in [T^. However, this factor affects the non-abelian, i.e. quantum 
gravity-dependent, structure of the theory in a very non-trivial way, and will affects the effective non-commutative 
field theory corresponding to it, as we are going to see in the next section. It is therefore important, and cannot 
be neglected. Of course, the main reason not to neglect it is that it follows directly and uniquely from the discrete 
path integral we have defined and computed, so it is a necessary part of the definition of a spin foam model for 3d 
quantum gravity. Also, let us stress that, had we used or obtained delta functions over the group enforcing the classical 
equations of motion, then this factor would have been fully equivalent to a constant normalization factor Cos (to), 
indeed present in the normalization used in and thus physically irrelevant. This is the effect of the ambiguity in 
the re-construction of the Feynman propagator from the Hadamard one that we had mentioned in section II. There 
is also a constant factor equal to g^^^^^-^ coming from every edge of F (also present in 0|). 

Putting everything together, we arrive thus at the following expression for the partition function for gravity coupled 
to matter 

J 11 9e' II 2 i-i- A{Sin{9,)+ie,f i-i~2Sin{m,) Sin^{9e)- Sin^{me)+iee 

e*6^A* e e^r v v / 

The properties of the above expression are very similar to those of (|^ . In particular 

• The amplitude is sharply peaked around the classical configuration, which is = for e ^ F, and 9,, = ifi^ for 
e G F, in agreement with one's expectations. 

• It is clear that H3U|) changes into the complex conjugate of itself, if one switches from integrating over the x's 
that make a positive dot product with appropriate vectors to those which a make a negative product with those 
same vectors. This is precisely what one would expect from continuum considerations. 



The model is not triangulation independent, essentially because the particle- free part of the amplitude is not 
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Since we chose the triangulation to be adapted to the particle graph, to start with, there is not much sense in asking whether the model 
is triangulation independent at the location of the particle (Feynman graph). 
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One can of course give a spin foam presentation of the above partition function. This is given by: 

^= E ll('^Je + l)l[K,Am,,e) n {6j} 

labellings e^F eGF tetrahedra^A 

The sum goes over all labellings of the faces of the dual two-complex by representations of SU(2). In the above 
{6j} denotes the usual 6-j symbol, and the Kj{eys are given by: 



2Sin{rhg) Si'n?{9) — Siri^^fhe) + it 



3 



2 Cns(O) 7 „• „-i(2i+l)()fi„-ie) 

F . \ t d. 

2nSin{rhe) J (a;^ + e2"")(a;2 + e-2"«) 

The Fj 's can in fact be evaluated explicitly in terms of elementary functions (which can be neatly repackaged into 
the hypergeometric functions), however what is important about them is that: 

• For half- integer representations (j is a half- integer), Fj is zero. This is immediate from the fact that the integrand 
is odd in this case. 

• Fj is manifestly real. 

Because of these two properties, the sum of Kj and its complex conjugate, i.e. the real part of the above 
coefficients, is simply given exactly by the character Xjl^™^ "'''), so that the sum of the whole partition function 
with its complex conjugate gives back exactly the usual coupled Ponzano-Regge model, as one expects. This has 
to be compared with the proposal for the Feynman propagator in pl| . where such correspondence is realized 
only modulo proportionality factors (functions of me) that are instead fixed uniquely from first principles in our 
procedure. 

It is very interesting to look at the quantum field theory (abelian) limit of H30|l . along the lines of We turn our 
attention to this task now. 



V. THE QFT LIMIT 

We want to see now how ()30|l compares to the evaluation of the usual Feynman graphs in fiat spacetime. This will 
imply first the re-writing of the above partition function, in the appropriate approximation, to be explained below, 
as the evaluation of the Feynman diagrams of a non-commutative field theory that takes into account the quantum 
gravity induced modifications on the propagation of matter only; second, it will imply checking that the same partition 
function reduces to the evaluation of the Feynman graphs of ordinary scalar field theory in fiat space in the no-gravity 
or G ^ limit. 

In other words, we want to concentrate on the matter part of H30() and 'integrate out' the gravity degrees of 
freedom. To do that note that by far the dominant contribution to the partition function from the no-particle 
edges, comes from the configurations which are fiat on those edges. This means that, as far as the particle edges are 
concerned, the only contributions that are relevant are those for which the connection is flat everywhere except at 
the particle graph. Of course, this is an approximation, we are replacing the fully quantum partition function, with 
the one in which gravity is fiat (classical, on-shell), everywhere where there is no matter. 

To turn the above into formulas, note that action (|10() can be written as 

S= f Tr{EAF)-y2f Tr{E p,) ^ Sp{g^,) + S,{g^,,p) 

JM j Jji 

where g^i, stands of course for the metric. We have switched to the second order formulation to make the physical 
point more transparent. Now, what we want to do is to approximate the partition function in the following way: 
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i i 

i i 

in the above gjjj, and g'^^ denote the metric that appears in the pure gravity or in the coupled part of the action 
respectively, rj^i, is the Minkowski metric.^" 

Putting the above differently, gravity has a dual role in this model. It has independent (only quantum) degrees 
of freedom itself, and it couples to and deforms (as we shall see below) the matter action. What we want is to 
freeze the (quantum) dynamics of gravity for now and concentrate on the deformation it causes in matter propagation. 

Keeping the above in mind, we replace the pure-gravity amplitudes in (|30|1 with delta functions, obtaining^"'^ 

•J e*GA* e^r eSF ^ ^ ' 

as we shall see below, H31|l can be thought of as a Feynman diagram evaluation of a noncommutative field theory. 
Since we will be interested in taking the limit in which we recover the usual, flat-space QFT, we must put Newton's 
constant G back into the equations. 

To see where the G should go recall that, as was mentioned in section II, fhe is the deficit angle associated to the 
mass me- More precisely, comparison with the classical solution reveals that 

Deficit Angle = ifie = AirGme — nnie 

Thus we have to replace everywhere mg with nnie- It is natural then to replace 0e everywhere with kO,,, and 
interpret the equation above as an equation between momenta^^ 

K Momentum flowing on edge e — k. Mass on edge e 

This reinterpretation is clearest on-shell, but we will retain this interpretation off-shell as well, thus interpreting the 
gravitational holonomy corresponding to 6 as the true quantum (and thus off-shell) momentum of the particle. As 
we have stressed above, this is of course consistent with what we know [T^llfl lT^llTj l about quantum point particles 
coupled to 3d quantum gravity, as well as with our own results. So, we shall parametrize our group elements as 

g^Cos{Ke)l + J-nSin{Ke) (32) 

The above parametrization is a very physically transparent one, and it is most useful for taking the k — > limit. To 
get the noncommutative held theory, we shall use however a slightly different way of representing the group elements 
given by 



g = K^P\n + kP ■ J. 

In doing so we have introduced the non-commutative momentum coordinates P, valued in B^, the 3-ball with radius 
1/k in M"^ , previously mentioned in section II and first used in . 



■^^ It may help to note that the discrete analogue s of Sp{g^v) and Sc{gij.iy,Pi) were given in I24i . 

This expression is exactly the one obtained in except for the Cosine factor. We follow the same normalization used there, dividing 
by the factor. 

This, in fact, is a highly nontrivial step which was first done in 14|, and although we are trying to make it look plausible, its real 
justification is an entirely a posteriori one: if one does this, then the usual quantum field theory is obtained in the k — > limit. 
Following the definition, the integration measures for g and P are related by: 



dg ■■ 



SU(2) JBi Ji_^2lp\2 
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Now, using the star product introduced in [T^, given by 

and the 'Fourier' transform naturaUy associated to the above product^** 

it can be shown, using the very same procedure used in 14], that H31|l can be rewritten as^^ 

/ n ^ / n ^2 ' s^nHJ.^^ Y~ n (★eea.e'^--^-) (33) 

The integrals over X„'s give delta functions, enforcing momentum conservation at the vertices. The expression that 
one gets is 

e •' eer Ke| ^ H ICe „gr 



,2\E\ 



It is clear that the above expression has exactly the form of a usual Feynman diagram evaluation, but in terms 
of the non-commutative momenta P, i.e. for a non-commutative field theory that we are going to write explicitly 
in the following. It looks the same as the one obtained in [T3. l2]| . except for the integration measure over the Ps, 
that is given here by an apparently usual k^(PP, while was given by dg — .'^ in p^l2ll|. It is also important 

^/ 1 — ^2 |p|2 

to notice that the fact that the (noncommutative) momentum is conserved follows here from moving from l|3()|l to 
(|31|l . In other words, conservation of momentum is only valid when gravity degrees of freedom are 'frozen', and the 
dynamics of pure quantum gravity as encoded in the partition function H30|l is neglected and the partition function 
itself is evaluated at the dominant contribution (flat spacetime) only. This means that momentum conservation, even 
at the non-commutative level, is only approximately satisfied in the causal theory. 

Inverting the expression for the Feynman propagator used in (|33|l . one can see that the action that generates 
Feynman diagrams of the form'^^ (|33|) is given by the following scalar noncommutative field theory action 

1 r (p^g)-^'^). ~ , X f 

S{<l>) = / dg , 4>{9)4){9 ) + ^ / dgidg2dg35{gig2g3)<j){gi)(t)ig2)(l}{g3) 



The coupling constant A has been introduced to make the action above look similar to the familiar 0"^ theory. 
The position space version of this action is given by: 



Sin (Km) _^ ^ A 



Vi - ^292 



0(a;)) + ^(^0*0*0) (x) 



(34) 



whereV is the Laplacian on M!^. Once more, the action differs, in both momentum and position space, from the 
one obtained in [3 due to the additional Cos{9) factor per edge that is present in our partition function, 
and necessarily so because it follows uniquely from the discretization of the BF partition function plus the causality 
restriction. 

We shall take now the k ~ G limit of H31|) . Using the first parametrization (|32|l , it is obvious that integration 
over the group SU(2) is replaced by integration over its Lie algebra su(2) ~ R^. 

dge. ^ d^pe* 

5(7(2) JK3 



X is an su(2) element, while gi,P2 are SU{2) group element s. A lso, the minus in the exponent is purely a matter of convention. The 
choice here is made so that the formulas agree with those in |14| . 
\E\ as before is the number of the edges of the triangulation. 
If one restricts to trivalent graphs. 
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Where we use the notation p = Oft. It is clear that with this definition, the group elements (|22() can be expanded in 
powers of k as 



Keeping the above in mind, and using the same techniques as in |l4j . in particular using the triangulation indepen- 
dence of the approximate particle-free part of the partition function, ensured by the use of delta functions over the 
holonomies of the connection, we see that the abelian limit of the full partition function H3U|) is given by: 



The equation above is almost exactly the same as the expression obtained in with the sole (welcome) difference 
in that the amplitudes on the particle edges are the Feynman propagators instead of the Hadamard ones, i.e. is 
equal to the evaluation of the usual Feynman diagram associated to the graph F in flat spacetime for a scalar field. 
We see that, as anticipated, the abelian limit of the causal coupled partition function is insensitive to the edge factor 
Cos{0) that marked the difference between the non-commutative Feynman propagator as uniquely obtained from the 
causal restriction of the BF partition function and the one used in . 

Let us also emphasize once more that the noncommutative analogue of (|35() , is not exactly the effective field theory 
coming from the full partition function. What it does correspond to is the theory in which instead of the (^sin{0}+ie)^ 
factors for the particle-free edges, one has delta functions. In other words, it is the theory in which gravity is on-shell 
away from the particles. It is clear that the off-shell partition function (|30|l is not exactly equal to the noncommutative 
field theory written above, because of the 'virtual' gravity degrees of freedom present^*. 



Let us summarize our results. We have shown what is the correct discrete counterpart of the causality condition that 
reduces BF theory to 3d quantum gravity in the Lagrangian setting, i.e. the restriction on the geometric variables 
that imposes on the one hand the correct dependence of the path integral on the spacetime orientation, and, on 
the other hand, realizes the correct 'off-shell' (with respect to the canonical Hamiltonian constraint) propagation of 
gravity degrees of freedom. We have constructed the corresponding causal spin foam model for pure 3d quantum 
gravity, and discussed its properties. We have then coupled matter to the model, along the lines of T^l, and derived, 
through a rigorous discretization and evaluation procedure, a causal coupled spin foam model for 3d gravity and 
scalar matter (the model can however be trivially extended to include spin degrees of freedom, whose coupling is 
not affected by the causality restriction). Then, we have shown how the causal model admits a physically motivated 
approximation corresponding to approximating the integral over pure gravity degrees of freedom with its dominant 
contribution, represented by the flat spacetime geometry. This approximated model has then been re-written as the 
Feynman evaluation of a Feynman graph for a non-commutative scalar quantum field theory, encoding fully the matter 
degrees of freedom plus the quantum gravity corrections to matter propagation, while the un-approximated model 
encodes also the quantum gravity corrections to flat space momentum conservation. The amplitudes so obtained 
involve non-commutative Feynman propagators, as one would expect. No ambiguity at all arises in re-writing the 
approximate partition function in terms of such Feynman evaluations, and everything follows uniquely from the original 
discretization of the causally restricted partition function. Finally, we have extracted, from the form of the non- 
commutative Feynman amplitudes, the unique non-commutative field theory action producing them in perturbative 
expansion, both in momentum and configuration space., and shown that both the action and the corresponding 
Feynman amplitudes reduce to the correct usual Feynman amplitudes for a scalar field theory in flat commutative 
space in the no-gravity limit. 

Let us now conclude by mentioning only some of the many possible uses and further developments of the results 
we have presented. 



At the level of the action, it is very easy to see that I34i reduces to the usual scalar field theory action in the limit k ^ 0, since the star 
product goes to pointwise multiplication in this limit, and the derivative factor in the denominator of the kinetic term drops out. 
This makes equation I3UI dependent on the triangulation chosen, and thus on extra information not contained in the Feynman graph F 
alone. This last point makes unlikely, in our opinion, although not ruled out that there exists a different non-commutative field theory 
that gives the full true partition function as the evaluation of its Feynman graphs. 





(35) 



VI. CONCLUSIONS AND OUTLOOK 
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We can summarize what we have achieved as the reahzation, in a discrete context and in a spin foam language, 
of the path integral or sum-over-histories idea for 3d quantum gravity coupled to matter fields, in the Lagrangian 
framework. One can then use this path integral for computing quantum gravity transition amplitudes and observables, 
being assured that the formalism correctly encodes causality and orientation dependence restrictions. However, it is 
maybe of more direct physical interest to compute the effects of quantum gravity corrections to matter scattering 
amplitudes, now that the unique and correct expression for the quantum gravity modified Feynman amplitudes for 
matter fields have been derived. 

There are several reasons, however, for believing that only a group field theory context would render the use of the 
above amplitudes truly physically meaningful (as far as 3d quantum gravity can be physically meaningful, of course). 
One obvious one is the triangulation dependence of the causal model, both for pure gravity and matter coupling. 
This triangulation dependence has to be expected, as we have discussed, since 3d gravity lacks local degrees of 
freedom only on-shell, i.e. only after the Hamiltonian constraint has been imposed, that is what causality prevents us 
from doing, in a sense. On the other hand, this triangulation dependence is still physically troubling, as it indicates 
that we cannot restrict ourselves to a fixed choice of simplicial complex on which to define our theory, as this implies 
necessarily a truncation of the degrees of freedom that are captured by our model. The natural way to overcome this 
limitation, that is also the only known way in the physical 4d case, is to pass to a group field theory formulation, and 
to generate our spin foam model within a perturbative expansion of a group field theory that has the interpretation 
of a quantum field theory of simplicial geometry 0, Q, I3 . On the one hand this would extend greatly the scope and 
context of our model, given that GFTs include in their perturbative expansion a sum over different topologies as well 
as different triangulations and that at the same time point towards possible non-perturbative results as well |3j . On 
the other hand, going to a GFT context would first of all solve the issue with triangulation dependence, as we said, 
but also provide our causal model with an even more natural role. As we have discussed briefly in the introduction, 
in fact, the need for an 'off-shell'propagation of quantum gravity degrees of freedom and thus for the un-symmetrized 
Lagrangian sum-over-histories, as opposed to the symmetrized one realizing just the projection onto physical states, 
becomes absolutely apparent in a 3rd quantized context, i.e. that of a true field theory formulation of quantum 
gravity. This is exactly what is provided by the GFT framework for simplicial quantum gravity and spin foam models. 

As mentioned, a proposal for a class of generalized group field theories, corresponding to causal spin foam models, 
in any spacetime dimension and signature, exists |24|. and can be specialized to the 3d case in a rather straightforward 
manner. The first thing to check, therefore, is whether our causal model can be seen as arising from this class of 
GFTs. If it does not, the modification of the generalized GFT formalism of |2^ that instead does the job should be 
sought for. Having at hand such a GFT model, the coupling of matter fields for it, along the lines of [13, has to be 
investigated. Only after this, we think, we can truly start tackling physical issues on the basis of the causal model 
we have constructed. 

Of course, true physics needs a 4-dimensional setting, and thus the main aim should be to work in the context 
of 4-dimensional quantum gravity. However, as it is well known, no consensus has been reached yet concerning the 
correctness of any specific 4-dimensional spin foam model, and the analysis of the proposed ones has still to progress 
past the initial stage. This is true also at the group field theory level. Concerning the issues that have interested us in 
the present work, i.e. the issue of causality in quantum gravity, and in gravity coupled to matter fields, this is not really 
dependent on the spacetime dimension. Therefore we expect our procedure and results to be of direct application 
and extension to any 4-dimensional spin foam model. In particular, if the 3d case can be handled successfully within 
the generalized GFT context of [2^, this would imply almost automatically that the same formalism is capable of 
implementing correctly causality restrictions in any dimension. Much more work is also needed to understand the 
coupling of matter fields in the 4d context; however, once again, our procedure and result would most likely be of 
direct application to any spin foam model of 4d quantum gravity coupled to matter that makes an explicit connection 
with classical continuum and discrete actions. An interesting testing ground would be, for example, the model for 
topological gravity coupled to string- and brane-like excitations investigated in [29|. 
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Appendix A 

Here we collect a few facts about the group SU(2) and its Lie algebra su (2) that we used in the text. 
The basis for the Lie algebra is Ji = zCTj, where Ui are the Pauli matrices. Thus 

= ( ) ' = ( * ) ' = J 

These obey 

\Ji 1 J j\ '^^ijk Jk 

Tr{J,) = 



From which it follows that 

J ■ a J ■ b = —a ■ b — J ■ {a X b) 
The last equation has been used in converting the actions we dealt with (Q, Hll|) ') to vector form. 
Using the exponential map, the general SU(2) group element can be written as 

e^" -^" = Cos{9)l + J- nSin{e) 
Where r? is a unit vector, and is a real parameter which goes from to tt. 

Appendix B 

Here we prove somewhat more rigorously equations p9l) and H2()f) . The proof and notation below are from [s^ . To 
begin with we define: 

• A cone in R'^, is a set C such that if x e C then Ax € C, VA > 0. We shall restrict our attention to acute cones. 
The simplest definition for an acute cone is that which states that a cone is acute if its intersection with its 
inverse through the origin is empty. 

• The cone conjugate to C, is the cone C* such that ^ e C* if ^ • y > , Vy G C. 

• The tubular region T*^ associated to the cone C is the subset of C^, which is given by 

^m? + iC ^ [z ^ X + iy: X eS? , y e C] 

We let C now be a very 'narrow' cone along the z-axis. In other words in polar coordinates 

C = (r, e, 4>) e with Q<e <5,r>Q 

where (5 is a small number (which wc will want to take to zero in the end, since what we want is the degenerate case 
in which C is just the z-axis). Note that C* is almost the entire half-space with z > in this case. We claim that 



Jc Jc'f 



where z £ T*^. The integral over n^^ can now be done giving (|19|l and 126() 



This n is diflferont from the one used in parametrizing the group elements in the main text. 
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The proof is very simple. It can be shown |33| that the left hand side of the equation above is a holomorphic 
function of z, and it is clear that the right hand side is holomorphic, because the denominator never vanishes'^". As 
is well-known, if two holomorphic functions agree on an interval, then they are equal everywhere. Thus we will be 
done if we can show that there's an interval in along which the two sides of the equation are equal, we choose the 
subset of given by the purely imaginary elements, i.e. we let z = iy , x = 0. Then 



cns^ Jo 



e-'-(y")r^drd'^7 



2i 



d n 

cns^ [y ■ Jo 
_d^n 

c 



e du 



q.e.d. 



Doing now the integral over n, one gets 



which coincides with H19|) and H2t)|l . 



2TTtan^( 5) 



The left hand side of the equation above is called the Cauchy kernel for T'-^ . It is important to realize that for us, 
only the values for z whose real part is along C are the ones that are relevant. This is because our prescription of 
integration is dependent on z. 



Appendix C 

Here we show the Riemann surface for the integrand in H28|l . and the contours leading to the Feynman propagator. 

The Riemann surface for the function that we have, is a two sheeted one (this comes from the square root, of 
course). There are two branch points, at 1 and at -1. To obtain the surface, imagine that one has two copies of the 
complex plane with a cut along the real axis joining the two branch points. In fact, to make things easier to visualize, 
imagine that we add the points at infinity, and consider two Riemann spheres with cuts. Then we should glue them 
in the usual criss-cross way typical of the square root Riemann surface. It is easy to see that the surface obtained is 
topologically a sphere'^-'^ . 

The Riemann surface for the — j- that we have is shown in the first figure. The meridian that is closer to the 

viewer, corresponds to the copy of [-1,1] which has negative square roots, while the one that is further away (behind 
the sphere) is the copy that has positive square roots. 

There are three distinct cases for our integral: 

A- 9 <m. 

B- m < 9 < TT — m. 

C- t: ~ m < 9. 

The contours are as shown in the figure^^: 



This is why we needed to 'open' the z-axis slightly. 

In fact there are punctures at the locations of the branch points, but these can be filled unambiguosly. 

Needless to say, all contours which can be deformed into each other without passing through singularities or branch points are equivalent. 
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FIG. 1: This is the Riemann surface associated with — 3-. The arrowed lines are the contours of integration in the three 

(1— u;2)3 

cases mentioned in the text. 

A- In the first case the contour just follows the 'negative' interval^^, i.e. both endpoints are chosen to lie on the 
copy of [-1,1] which has negative square roots. 

B- In this case the contour follows the 'positive' interval from —Cos(6 — m) to zero, then moves to the 'other' zero, 
then continues along the 'negative' interval to —Cos{6 + m). 

C- In the last case, the contoiir follows the 'positive' interval. 

A somewhat more conventional picture of the contours is shown in the second figure. Here, the thick line joining 
the branch points is the cut, and as one approaches it from above one gets positive square roots, while approaching 
it from below gives the negative ones. 



One could reverse the conventions, by letting the first case be along the positive interval, and vice versa for the other two cases. This 
however introduces an overall minus sign into the final expression. We chose this prescription to maJje the abelian limit agree with usual 
convention for Feynman propagators. 
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Imaginary Axis 
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FIG. 2: Another picture of the contours of integration. As one approaches the cut (thick line joining -1 and 1) from above one 
gets positive roots, while approaching from below gives negative ones. 
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